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ABSTRACT 

The properties of elliptical galaxies are broadly consistent with simulated remnants of gas-rich 
mergers between spirals, motivating more detailed studies of the imprint of this formation mechanism 
on the remnant distribution function. Gas has a strong impact on the non-Gaussian shapes of the 
line-of-sight velocity distributions (LOSVDs) of the merger remnant, owing to the embedded disk 
that forms out of the gas that retains its angular momentum during the merger, and the strong 
central mass concentration from the gas that falls to the center. The deviations from Gaussianity are 
effectively parametrized by the Gauss- Hermite moments /13 and h±, which are related to the skewness 
and kurtosis of the LOSVDs. We quantify the dependence of the (/13, hi) — v/a relations on the initial 
gas fraction f gas of the progenitor disks in 1:1 mergers, using Gadget-2 simulations including star 
formation, radiative cooling, and feedback from supernovae and AGN. For f gas < 15% the overall 
correlation between ft, 3 and v/a is weak, consisting of a flat negatively correlated component arising 
from edge-on viewing angles plus a steep positively correlated part from more face-on projections. The 
spread in v/a values decreases toward high positive /14, and there is a trend toward lower /14 values 
as fg a s increases from to 15%. For f gas > 20% the (ft.3, /14) — v/a distributions look quite different 
- there is a tight negative h^ — v/a correlation, and a wide spread in v/a values at all hi, in much 
better agreement with observations. Re-mergers of the high-/ sas remnants (representing dry mergers) 
produce slowly rotating systems with near-Gaussian LOSVDs. We explain all of these trends in terms 
of the underlying orbit structure of the remnants, as molded by their dissipative formation histories. 
Subject headings: methods: n-body simulations - galaxies: elliptical and lenticular, cD - galaxies: 
formation - galaxies: interactions - galaxies: kinematics and dynamics 



1. INTRODUCTION 

Many of the observed properties of elliptical galax- 
ies indicate a violent formation history. They are 
dynamically hot systems, with high velocity disper- 
sions dominating over ordered stellar streaming. Gas- 
rich tidal tails, and rings and shells indicative of the 
recent disruption of a spiral galaxy, often surround 
systems otherwise resembling o rdinary giant ellipticals 
(|Arp| 119661: Ivan Dokkuml 120051). These observations 
first led " iToomre fc Toomrd l| 19721) to suggest that el- 
liptical galaxies are the products of mergers between 
spirals. N-body simulations of mergers between disk 
galaxies generically reproduce many of the gross prop- 
erties of observed ellipticals, including density profiles 
that follow the V 1 / 4 " law in projection (Ivan Albadal 
198llMcGlvnnlll98l. fl at rotation curves (|Whitdll978t 
Farouki fc Shapiro! I1982D. slow rotation and aniso tropic 
velocity distribu tions dWhitd [l978t iGerhardl 119811) . and 
triaxial shapes ([Gerhard! 119811 : iBar ncs 1988). Today 
we understand the Toomre merger hypothesis within 
the broader context of hierarchical structure forma- 
tion in the Universe - in the favored ACDM cosmol- 
ogy, an elliptical galaxy-sized halo has typically suf- 
fered ~1 major merger since z ~ 2 — 3, the epoch 
over which the bul k of its stellar mass forme d (e.g . 
LacevfcCoiel 1x993 IManriaue fc Salvador-Sold HoM 
Kauff mannl 119961 : iDe Lucia et all 120061 : iHopkins et all 



I2008ai l2009al IStewart et all 120081 ) . 

Observations of elliptical galaxies also show signa- 
tures of dissipation in their buildup. The phase space 
densities in the cores of ellipticals exceed the maxi- 
mum densities in spirals, implying a violation of Li- 
ouville's theorem if ellipti cals arise from dissipation- 
less spiral-spiral mergers (|Carlberdfl986T) . This prob- 
lem is naturally resolved by a gas component that is 
driven to the center by tida l torques during the merger 
(iNegroponte fc Whitd Il983t iBarnes fc HernquistJ 119911 : 
iHernquist et al.lll993l ). Elliptical galaxies also display 
kinematic features that are most easily explained as 
signat ures of dissipation, such as rapid oblat e rotation 
(e.g. IDavies et all 119831: IBender et all 1 19921) . embed- 
ded disks (IBenderi 119901: iRix fc Whitd 119901). and kine 
matically decoupled subsystems (IHernquist fc Barnei 
19911: lEfstathiou et al.lll982t iFranx fc Illingworthl 1198 
Jedrzeiewski fc_Schechterl 119881 : iKrainovic et alj 120081) 



Cox et alT i 20061 ) showed that a wide variety of the pho- 
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tometric and kinematic properties of observed ellipti- 
cals, including their half-light radii and velocity disper- 
sions, flattening, isophotal shapes, velocity anisotropy, 
and major and minor axis rotation were much better 
reproduced by an ensemble of simulated 40% gas disk 
merger remnants than by the corresponding ense mble of 
dissipationless mergers. iRobertson et al.l (|2006| ) found 
that a gas fraction of ~30% was needed to match the 
observed tilt in the "fundamen tal plane" ( [Dressier et al.l 
H987tlDiorgovski fc Davislll987D . a value close to the typ- 
ical gas fraction at z ~ 2 — 3. 

The violent relaxation in galaxy mergers is incomplete 
- substantial memory of the initial conditions is retained. 
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For example the remnant may tend to rotate in the sense 
of the initial orbital angular mome ntum or disk spins 
(|Barneslll992l: lHernquistlll992l 119931) . and an initial stel- 
lar metallicity gra dient will be blurred but not erased 
during the merger (White 1978). This incompleteness is 
enhanced by the presence of gas, which can form cold 
features in the stellar distribution relatively late in the 
merger process, which subsequently experience less vio- 
lent relaxation than features present at the beginning of 
the merger. This fine structure provi des a fossil re cord 
of the galaxy's formation history (e.g. IStatlerlll987l ). 

One way to parametrize the finer features in a galaxy's 
distribution function, which is well-suited to spectro- 
scopic observations, is through the moments of the line- 
of-sight velocity distribution (LOSVD) as a function of 
projected location in the galaxy. For continuous stellar 
line profiles, these moments are measured by fitting the 
LOSVD to a Gaussian multiplied by a truncated series 
of Gauss-Hermite (GH) basis functions Hk(w), 

n 

F GH (v los ) X e-^ w2 [l + J2hkH k (w)} : (1) 

where w = (vi os — v)/a and a is the line-of- 
sight (LOS) velocity dispersion. This method is 
less sensi tive to noise than computing the moments 
directly (Ivan der Marel fc Franxl [1991 : iGerhardl H991 : 
IBinnev fc Merrifieldl 119981 ). For modest deviations from 
Gaussian ity, h% represents the skew ness of the distribu- 
tion (see IBinnev fe Mer rifield 19981 Fig. 11.5) and can, 
for example, indicate whether the net rotation in a galaxy 
arises from the dominant orbit population, or a small 
group of streaming stars. The quantity hi measures the 
kurtosis; positive hi values generally indicate a radial 
anisotropy in the velocity dis tribution, whil e negative hi 
indicates a tangential bias (jGerhardl Il993f ). In a non- 
rotating system, hi can distinguish between a dominant 
population of radial (box) orbits, and cancelling streams 
of high-angular momentum (tube) orbits. 

In recent years, integral field spectroscopy (IFS) has 
made it possible to obtain high S/N, high-r esolution 2D 
maps of the LOSVD i n nea rby ellipticals ([Bacon et al.l 
120011: lEmsellem et al.l I2004D . Schwarzschild modelling 
studies have shown that, in practice, the observed 2D 
maps of the first four moments of the velocity distri- 
bution, {v, a, /13, hi}, are sufficient to uniquely recon- 
struct the full 3D stellar orb i tal distribution in most 
cases ([CappeHari et al.1 120031: Ivan de Ven et al.l 120081: 
van den Bosch et al]|2008l : Ivan den Bosch fc van de Venl 
20081) . Patterns in the higher moments, {ht,k > 5}, 
typically fall below the noise in the observations. It is 
therefore standard practice to truncate the GH series at 
k = 4, and encapsulate the non-Gaussian shape of the 
distribution in the t wo parameters and hj. 

Previous authors (iBendo fe Barnes! 2000: J esseit et al.l 
120051: INaab et alJ2 006t l Jesseit et~a!l 2007) have observed 
that simulated gas-rich merger remnants occupy different 
areas of the ^.3,4 — v/a planes from dissipationlcss disk- 
disk merger remnants. Dissipationless remnants tend 
to show an overall positive correlation between /13 and 
v/a, while in gas-rich remnants, as in observ ed ellipti- 
cals ([Bender et al.lll994l:lKrainovic et al1l2008f). th is cor- 
relation is negative. iGonzalez-Garcia et al.l (|2006f ) show 
that adding a bulge component to the merging disks can 



also yield larger /13 values and a negative /13 — v/a cor- 
relation, by allowing the disks to retain more of their 
initial angular momentum. I Bender et alj (j!994f ) demon- 
strate that the observed h% — v/a relation is too steep 
to be explained by two-integra l oblate rotator models 
(e.g. iDehnen fe Gerhardl fl994) . but can easily be ac- 
counted for by the superposition of a hot stellar spheroid 
with a cold embedded disk. This type of distribution 
function arises naturally from a combination of violent 
relaxatio n and dissipation in gas-rich mergers of spiral 
g alaxies dHopkins et al.ll2009bD. 

Uesseit et al.l ([20051 ) and INaab et all (j2006f ) performed 
a detailed study of the stellar orbit structure of 1:1 and 
3:1 merger remnants with and 10% gas, and its imprint 
on their photometric and kinematic properties. They 
found that the gas drives an exchange between box and 
short-axis tube orbits, making the remnants more oblate 
in shape. This suppression of box orbits strongly influ- 
ences the shapes of the LOSVDs, bringing the h 3 ^ — v/a 
relations into bett er agreement with ob served rapidly ro- 
tating ellipticals. J esseit et al.l ([2007D showed that 2D 
kinematic maps of the same set of remnants display many 
of the intriguing features seen in real galaxies, including 
misaligned rotation, embedded disks, and kinematically 
decoupled cores. 

In this paper we explore the effect of gas on the 
^3.4 — v/a relations using a version of the tree/SPH code 
Gadget-2 (Springcl 120051 ) that includes star formation, 
radiati ve cooling, and feedback from supernovae an d 
AGN (ISoringel fc HernmiisH [20031 ISnringel et al.l 12005?) . 
The inclusion of real-time star formation allows us to con- 
sider the higher gas fractions typical of spiral galaxies at 
the peak elliptical formation epoch, f gas ~ 30% at z ~ 2. 
For fixed f gas , the dissipative signature on the dynamics 
may be reduced by star formation, if the gas is converted 
to collisionless material early on in the merger. The dis- 
sipational features may also be more spatially extended 
when star formation is included. 

We simulate a representative series of mergers between 
equal-mass, Milky Way-sized disk galaxies at a series of 
initial gas fractions ranging from to 40%, as well as 
re-mergers of the spheroidal remnants. We quantify the 
dependence of the /13 4 — v/a diagrams on f gas , and ex- 
plain the patterns that we find in terms of the underlying 
stellar orbital distribution. Our primary goals are: (a) 
to see whether this diagnostic points to the same typi- 
cal gas fraction as other indicators, such as the FP tilt 
and peak quasar redshift; (b) to locate "wet" (gas-rich 
disk-disk) and "dry" (re-mergers of gas-poor spheroid) 
mergers in (/13.4 — v/a) space, with the aim of distin- 
guishing these two populations in IFS observations; and 
(c) to connect the observed trends in the non-Gaussian 
moments with the underlying orbit structure typical of 
merger remnants in an intuitive way. 

2. SIMULATIONS AND METHODS 

Our initial galaxy models consist of exponential disks 
embedded in Hcrnquist (1990) dark matter halos, with 
masses comparable to th e Milky Way. The gal axy models 
are described in detail in lSpringel et al.l ((2005) . The halo 
concentration and spin parameter were set to c = 0.9 and 
A = 0.033. The disk mass was 4.1% of the halo mass, 
and its specific angular momentum was assumed equal 
to that of the halo. The halo was realized with 120,000 



3 



particles, and 80,000 stars initially comprised the disk. 
A fraction f gas of the disk stars were replaced with SPH 
gas particles at the start of the simulation, and this com- 
ponent could subsequently form new stars if a threshold 
density p t f,, tuned to match the Schmidt Law ((Schmidt 
1959: IKennicuttl [T998). was exceeded. The gravitational 
softening length was 140 pc, which set a minimum spa- 
tial resolution scale in the core. (See also lCox et aHl2006l 
for further description of the simulations.) 

The disks began at a separation of lOOkpc, on parabolic 
orbits with a periapsis distance of 7.1kpc. They typically 
merged after around 1.5Gyrs, and the simula tions were 
run to a final time t en d = 3/i -1 =4.3 Gyrs. ICox et aTl 
( 2006) show that the global remnant properties reach a 
steady state at one effective radius (R e ) around 0.3Gyrs 
after the merger, so by t en d the remnants can safely be 
considered relaxed at this scale. The set of eight merger 
trajectories used in this study correspond to orbits i— p in 
I Cox et all (|2006l ). chose n to uniformly sample the phase 
space of possible orbits (|Barneslll992( ). 

The scries of eight representative disk-disk mergers was 
performed at seven different gas fractions, f gas =0, 5, 
10, 15, 20, 30, and 40%. We also ran two series of re- 
mergers of the 20 and 40% gas remnants, intended to 
represent dry mergers of realistic gas-poor ellipticals. To 
produce a representative sample for each re-merger se- 
ries, we randomly selected eight pairs from among the 
disk-disk merger remnants, and re-merged them on tra- 
jectories i — p again. 

For each of the simulated remnants we selected 100 
random isotropically sampled viewing angles, and a syn- 
thetic 2D map of the kinematics was constructed within 
lR e for each LOS. A constant mass-to-light ratio was 
assumed throughout the analysis. In each spatial bin 
we constructed a velocity histogram and extracted v, a, 
and /14 by performing a least-squares fit to the 5- 
parameter function 



F{y) = Aexp[- 



2a 2 



a 



using the Levenberg-Marquardt method (|Press et al.l 
11992ft . We used 40 x 40 spatial bins and 80 veloc- 
ity bins within ±3.5w rms . The Gadget particles were 
smoothed over a radius h smoot h — max(h see ,h ng b), where 
h see = 150pc corresponds to a seeing of 1.5" at 20Mpc, 
and h ng b is a smoothing length based on the distance to 
the 128th nearest neighbor. This ensured that the effec- 
tive particle count in each spatial bin was > 1000. At t en d 
we also froze the potent i al, ex panded it in the basis set 
of iHernquist fc Ostrikerl (|1992f L and classified the orbits 
of the stellar particles in this potential. The results of 
this e xercise are thoroughly presented in (Hoffman et al. 
l2009h . 

3. RESULTS 

3.1. The h 3 — v/cr distribution 

Our results for the h 3 — v/a relation are presented in 
Fig. 1. The shape of the distribution is quite sensitive 
to fgas for the disk-disk merger remnants. At low gas 
fractions the overall correlation between h 3 and v/cr is 
weak; the distribution has a diamond shape consisting of 



a flat, negatively correlated component plus a steep, pos- 
itively correlated one. Although the intrinsic structure 
changes rapidly from flattened, prolate-triaxial systems 
dominated by box orbits to rounder, more oblate systems 
dominated by short axis (z—) tube orbits as the g as frac- 
tion increases from to 15% ((Hoffman et al.ll2009f ). there 
is no strong signature of this transformation in the asym- 
metric h 3 moment. 

For fg as > 20% the shape of the h 3 — v/a distri- 
bution changes rather abruptly. These gas-rich rem- 
nants are characterized by high h 3 moments that are 
strongly anti-correlated with v/a, i n better agreement 
with observed fast-rota tor ellipticals (|Bender et al.lH994l : 
iKrainovic et al.l l2008f) . The observations thus favor 
typical progenitor gas fractions > 25%, in agreement 
with other indicators such as the sha pes o f the rem- 
nants (|Cox et al.l l2006b iHoffman etHI 120091 . the two- 



component stru cture of ongoing merg ers ([Hopkins et al.l 
l2008bD. cuspv (IHopkins et all l2009d ) and cored ellip- 
ticals dHopkins et alj l2009dl), and the tilt o f the FP 
([Robertson et all 120061 : IHopkins et alj l2008cft . In the 
SAURON data there is also evidence for a second popula- 
tion of fast rotators with low h 3 values (|Krainovic et all 
120081 ): these systems cannot be explained by the simula- 
tions performed in this study. 

Re-mergers of the high-gas fraction remnants still show 
a distinct negative correlation between h 3 and v/a, but 
the remnant LOSVDs are much more concentrated near 
the origin of the h 3 — v/a plane. Note the much lower 
spread in re-merger h 3 values at v — compared to disk- 
disk mergers of any gas fraction. This result is not sur- 
prising, since violent relaxation in the re-merger tends to 
blur out structure in the distribution function and drive 
the LOSVDs closer to a thermal distribution. Fig. 1 
quantifies the degree of thermalization in a single dry 
merger. 

In Fig. 2 we illustrate how the positive and negatively 
correlated components arise in the low-/ gas remnants, 
and how the strong h 3 — v/a anti-correlation arises in 
the gas-rich ones. In dissipationless remnant I (top six 
panels), the major-axis (a;—) tube orbits produce sub- 
stantial minor axis rotation even though they only com- 
prise ~20% of the stellar mass within lR e , since they 
are highly streaming. In a pixel near the minor axis in 
an edge-on projection (looking along the intermediate 
axis), x— tubes dominate the luminosity, and when com- 
bined with the z— tube and box orbits peaked at v = 
they result in a distribution with an extended tail in the 
direction opposite the mean velocity. Since the veloc- 
ity distribution of the :r-tube orbits is not too strongly 
peaked, the /13 value is not very high and the h 3 — v/a 
relation is relatively flat. 

When the same remnant is viewed face-on (projected 
along the short axis), the three orbit classes are all piled 
on top of each other, and the streaming x— tubes form a 
bump in the tail of the dominant distribution of boxes 
and z— tubes. We thus get an extended tail in the same 
direction as the mean velocity, i.e. a positive h 3 — v/a 
correlation. This high-w bump in the tail produces a 
stronger asymmetry in the LOSVD while the velocities 
in the face-on projection tend to be lower, so this effect 
produces a much steeper h^ — v/a relation than the effect 
illustrated in the top row. 

The locus of face-on and edge-on projections in the 
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Fig. 1. — Local /13 — v/cr relation for each of our nine merger simulation sets. The first seven panels are for disk-disk mergers of different 
gas fractions, while the last two panels (lower right) are for dry mergers of the 20 and 40% gas remnants. Each map is averaged over the 
eight unbiased merger orbits, and f 00 isotropically selected lines of sight per remnant. The shading is proportional to the logarithm of the 
luminosity in each /13 — v/a bin, normalized so that the entire map sums to one. The red lines are contours containing 68 and 95% of the 
luminosity. 



hs—v/cr diagram, averaged over all remnants of a given 
gas fraction, is shown in the bottom row of Fig. 2. To 
illustrate the trends we show one set of low-/ 5as disk-disk 
mergers, one set of high-/ ffas disk-disk mergers, and one 
set of re-merger remnants. As indicated in the above ex- 
ample, face-on projections dominate a steep positively- 
correlated component, while edge-on projections domi- 
nate the shallower negatively correlated component. The 
positively correlated component is relatively suppressed 
at high gas fractions. 

In gas-rich remnants the dominant population (z-tube 
orbits) is streaming as in the edge-on dissipationless rem- 



nant, only now the streaming is much stronger. This ef- 
fect arises from the embedded disk that re-forms after 
the merger out of the gas that d id not lose its angula r 
momentum and fall to the center ([Hopkins et al. 2009b). 
The red and blue dashed lines in the bottom LOSVD 
show the contribution from "old" stars put in as colli- 
sionless particles at the beginning of the simulation, and 
from the 15% of the new stars (formed out of gas during 
the merger) that formed latest. These late-forming stars 
comprise only a few percent of the total stellar mass, but 
are prominent in spatial bins lying along the major axis. 
These disk stars have a much lower velocity dispersion 
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Fig. 2. — Explaining the /13 — t;/cr trends. Top row: Merger orbit I, 0% gas, at a nearly edge-on viewing angle. Second row: Same 
remnant, viewed nearly face-on. Third row: Merger orbit m, 40% gas, viewed at an angle about 30° from edge-on. First two columns: v 
and /13 maps within lR e . The color scale runs from —h max to h max , where h max is the maximum value of the moment in any pixel, and 
the number in the upper right corner gives the value of hmax- The black contours are over-plotted isophotes. The third column shows the 
smoothed LOSVD at the location specified by the white dot on the v and /13 maps. The heavy black line is the full LOSVD, while the 
red and green circles and black squares break it down by orbit class. In the case of the gas-rich remnant (bottom row), we also show the 
separate contributions from old stars (red dashed line) and the 15% of the new stars that formed latest (blue dashed line). (See the text 
for further discussion of these plots.) Bottom row: Net contribution of face-on and edge-on projections to the /13 — v/cr distributions. The 
three maps correspond to the second, sixth, and ninth panels of Fig. 1, but now each bin is colored according to the luminosity-weighted 
mean | cos0|, where 8 is the polar angle between the LOS and the intrinsic short axis of the triaxial remnant. | cos0j = 1 means perfectly 
face-on, while | cos#| = is perfectly edge-on. 
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than the remnant as a whole, and produce a sharp high-t; 
peak in the dominant z— tube population, giving a dis- 
tribution that is strongly skewed in the sense opposite 
the mean velocity. 

As a rule, disk-disk mergers tend to produce large 
asymmetric deviations from Gaussianity owing to sub- 
components that retain memory of the initial streaming 
of the disk stars, and to gas that is redistributed and spun 
up through dissipation to form a thin embedded disk. 
The /13 — v/a correlation is generally negative when the 
dominant orbit population tends to stream, and positive 
when the rotation is caused by sub-dominant streaming 
population superimposed on the tail of the dominant or- 
bits. 

3.2. The hi — v/a relation 

We now turn to the symmetric deviations of the 
LOSVDs from a Gaussian shape, measured by the "kur- 
tosis" parameter /14. In Fig. 3 we present the /14 — v/a 
relations from our simulations. For f gas < 15% the 
/14 — v/a distribution has a triangular shape that nar- 
rows sharply toward high positive /14. Note that this 
is the naive expectation for a single population of tube 
orbits - viewed edge on the high-angular momentum or- 
bits tend to pile up in the tails of the velocity distri- 
bution and any significant imbalance in the sense of 
their rotation produces a high v/a, while viewed face- 
on they pile up in a peak around v = 0. However 
this pattern differs sharply from IFS observations, in 
which a mild inc rease in the v/a spread is actually seen 
toward high /14 (jKrainovic et al.|[2008l ). There is also 
a trend toward lower hi values as f gas increases from 
to 15%, arising from the conversion of box orbits 
into z— tubes by the central mass concentration formed 
by the g as that falls to the center and produces the 
starburst (iBarnes fe Hernquistl [l996: Mih os fc Hernquisll 
fl99ilNaab et al.H2006tt . " 

We get a marked change in the shape of the hi — v/a 
distribution at f gas ~ 20%, as we did for the /13 — v/a 
relation. At higher gas fractions we get strong positive 
as well as moderate negative /14 values over a wide range 
in v/a - the v /a spread actually increases mildly toward 
high hi as in observed systems. 

Like the transformation in the /13 — v/a distribution, 
this rapid transition owes to the appearance of a strong 
disk component starting at f gas ~ 20%, as shown in Fig. 
4. In high-v/er, edge-on projections of the disk, /14 is in- 
fluenced by two competing effects - the strong negative 
kurtosis of the more isotropic population of old stars on 
z— tube orbits, and the strongly peaked population of 
late-forming disk stars that tend to produce a positive 
/14 in superposition with the old component. In the pixel 
shown in the top row of Fig. 4, the latter effect wins 
out. Although the old stellar z— tube population has a 
negative kurtosis, the narrow superimposed peak from 
the late-formed disk gives the full LOSVD a net positive 
/14. This positive /14 is greatly enhanced for the z-tube 
orbits alone. In the outer dotted pixel whose LOSVD is 
shown in the middle row, the disk component is some- 
what weaker and the negative kurtosis of the old z-tube 
population wins, giving the full distribution a net nega- 
tive /14. 

This example demonstrates how we may get a wide 
spread in hi , including strong positive as well as negative 



values, at high v/a. We caution the reader that in this 
situation hi is quite sensitive to small changes in the 
disk strength, and the GH fitting is sensitive to noise, so 
care must be taken in quantitative comparisons between 
simulations and observations. 

In the hi — v /a diagram the re-mergers are once again 
more concentrated around the origin, for the same reason 
discussed in the previous subsection. The smaller non- 
Gaussian moments of the re-mergers at fixed v /a suggest 
one way to separate re-merger remnants from other slow 
rotators. The re-mergers tend to have LOSVDs with pos- 
itive hi, also an expected signature of violent relaxation, 
which ge nerally results in a radiall y-biased distribution 
function (Bin nev fe Tremainell2008T ). 

Face-on projections generically give high positive hi 
values, as both the box and z— tube orbits are strongly 
peaked around v — in this projection. Particularly 
strong positive hi values can arise in face-on views of gas- 
rich remnants since the latest-formed stars form a narrow 
peak centered on v — 0, as illustrated in the third row of 
Fig. 4. The bottom row of this figure shows the average 
locus of edge-on and face-projections in the v/a plane - 
face-on projections dominate the high-/i4, low v/a part 
of the diagram, while edge-on projections dominate the 
high v/a and low /14 regions, as well as the high v/a, 
high hi wings in the gas-rich remnants. 

4. DISCUSSION AND CONCLUSIONS 

Mergers between spiral galaxies generically produce 
remnants with significantly non-Gaussian LOSVDs, ow- 
ing to the presence of subcomponents that retain mem- 
ory of the initial conditions. This effect is especially pro- 
nounced in mergers with f gas > 20% where the dissi- 
pative component re-forms an embedded disk, produc- 
ing a strong negative h^ — v/a correlation and a wide 
spread of hi values at high v/a. Both of these features 
are in agreement w i th observations of elli ptical galaxies 
(jBender et al.l H99l [Krai novic et al.|[200l . 

Upon closer examination the SAURON fast-rotator 
population appears to consist of two subgroups, one with 
a steep h^ — v/a correlation, an d a second with low /13 for 
all v/a (jKrainovic et al.H2008f ). A gas- rich major merger 
scenario cannot account for the second group, which sug- 
gests that the two fast-rotator populations may differ in 
their formation mechanisms. 

Dry mergers produce slowly rotating remnants with 
small /13 and /14, as violent relaxation blurs features in 
the distribution function and partially thermalizes the 
local velocity distribution. This provides one way to dis- 
tinguish dry merger remnants within the population of 
slowly rotating ellipticals. The large s pread in h% and 
hj am ong the SAURON slow rotators (jKrainovic et al.l 
2008) suggests that the slow rotators are not all dry 
merger remnants, as has been proposed to explain 
the fast-slow rotator dichot omy (jGonzalez-Garcia et al.l 
[200llKormendv et al.ll200l . 

At fixed redshift a tight correlation is observed between 
the gas fraction and stellar mass of late- type galaxies, 
with f gas varying from ~50% at M star = 10 10 M Q to 
~10% at M star = 3x 10 11 M Q at z ~ 2 (| Stewart et all 
120091 and references therein). The f gas labels in Figs. 1 
and 3 may therefore be taken loosely as a proxy for galaxy 
mass, and our gas-rich merger model predicts a variation 
in the /13.4 — v/a relations with galaxy luminosity. 
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Fig. 3. — Local — v/a relation for each of our simulation sets. The shading and contours are as in Fig. 1. 



In this paper we considered only a very restricted set 
of merger models, namely binary 1:1 mergers of bulgeless 
disks. The embedded disk component in our remnants 
may be artificially pronounced because real galaxies 
would suffer subsequent harassment through minor merg- 
ers, interactions with satellites, and secular processes. 
Some of the effects of increasing f gas (e.g. higher v/a 
and /13) ca n also be achieved by decreasing the merger 
mass ratio dBournaud et al.ll2005tlGonzalez-Garcia et al.1 
[20061: iNaab et all 120061 : Uesseit et all 120071 120091 ). and it 
is important to disentangle these trends. Extensions of 
this work include varying the mass ratios and bulge frac- 
tions, simulating cosmologically-motivated merger se- 



quences (|Burkert et al.ll2008l : IVailuri et al.ll2009f h and in- 
cluding secular heating processes. 

Ultimately we hope to make rigorous statistical com- 
parisons between large libraries of simulations and IFS 
data. This will require a careful consideration of obser- 
vational selection biases as well as systematics associated 
with the GH fitting procedure in the presence of obser- 
vational noise as well as particle discreteness in the sim- 
ulations. This study motivates such work by illustrating 
how detailed 2D kinematic information can reveal the un- 
derlying orbit structure of elliptical galaxies, and place 
stronger constraints on models of their formation. 
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/14 maps for this remnant. Second row: /14 maps for the z-tube orbits alone. The left map is for all of the stars on z-tube orbits, while the 
right map includes only old stars on z-tube orbits, excluding those that formed from gas during the merger. In the right-hand column the 
LOSVDs are shown at the two points indicated by the white dots, broken down by orbit class and stellar age as in Fig. 2. Third row: /14 
maps for merger orbit p, viewed face-on, at 0% gas (left) and 40% gas (right). The LOSVD at the dotted location on the 40% gas map is 
shown in the third column. Bottom row: Net contribution of face-on and edge-on projections to the /14 — vja distribution, plotted as in 
the bottom row of Fig. 2. 
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